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Abstract 

Using effective- field theory with correlations, we investigate the effects of interfacial pseudo-spin 
coupling fluctuations on the susceptibility and polarization of ferroelectric superlattices within 
the framework of transverse Ising model. It is found that the interfacial coupling fluctuations 
increase the susceptibility in the low temperature region. For a strong interfacial coupling, the 
phase transition temperature decreases with the strength of fluctuations of the interfacial coupling. 
The dependence of the susceptibility on the superlattice period of BaTiO^/ SrTiOs are plotted 
for different interfacial coupling fluctuations strength. At room temperature, when the interfacial 
coupling fluctuation increases, the peak position of the susceptibility will shift to a large superlattice 
period. 
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I. INTRODUCTION 



Recently developments in the fabrication of thin films have been applied to the ferro- 
electric thin films. Due to their good ferroelectric, dielectric, piezoelectric and pyroelectric 
properties, ferroelectric heterostructures have attracted much attention and have many po- 
tential applications, for instance, to nonvolatile dynamic random access memories(FDRAM), 
thin film capacitors, detectors, sensors, optical instruments. New ferroelectric materials with 
excellent dielectric properties at small sizes have also been considered to make electric device 
such as the small size capacitors. Thus, a high dielectric constant film with thickness less 
than O.lyum is a target for ferroelectric research. 

Giant permittivity associated with the motion of domain walls was reported in epitaxial 
heterostructures having alternating layers of ferroelectric and nonferroelectric oxides Q|. 
Experiments have also been carried out on the dielectric enhancement and Maxwell- Wagner 
effects in ferroelectric superlattice structure [^. The lattice mismatch and interfacial strain 
in ferroelectric multilayer are thought to be the main causes for the dielectric enhancement 

m 

. Theoretically, it is found that the dielectric property of the ferroelectric superlattice is 
very sensitive to the interfacial coupling and the thickness of the component 0] . The effects 
of stress on ferroelectric thin films have been studied within the framework of Landau theory 
m .he co„cl.,o„ ... M,.e. ten.,e ... enhance. wMle the M.he. co,„p.e.ve 
stress reduces the mean susceptibility. By taking the four-spin interaction into account, a 
: 'erroelectric-ferroelectric phase transition was found for large periods at low temperature 
3]. In Refs. y, 1^, the effects of long-range interactions and a non-ferroelectric layer on 
the dielectric properties of ferroelectric multilayer were studied. Defects zones and the 
mismatch at the interface between successive ferroelectric layers were observed in Ref. P|. 
It is anticipated that a ferroelectric multilayer with a high concentration of interfaces (where 
the bonding and the structure will in general depart from that of the interior of the layers 
or the bulk)will have some new properties. A complex and inhomogeneous interface may be 
important, or even dominant, in the superlattice structure when the individual constituent 
layer is only a few unit cells. So far, the effect of the interfacial structure on the dielectric 
property of a ferroelectric superlattice has not been thoroughly investigated. Analogous to 
the structure fluctuation or the bond randomness in the amorphous ferromagnets [10], the 
disordered structure of the ferroelectric interface must give rise to the randomness in the 
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interaction of the dipolar moments. 

In this paper, we introduce interfacial pseudo-spin couphng fluctuations in the transverse 
Ising model (TIM) to investigate the effects of interfacial structure inhomogeneity on the 
dielectric properties of ferroelectric superlattices. The coupling fluctuations may also exist 
in the interior layers of the superlattice. However, compared with the magnitude of the 
coupling fluctuation within the interface, the fluctuations in the interior layers are weaker. 
We focus our attention on the effects of the interfacial pseudo-spin coupling fluctuation. The 
dielectric property and phase transition temperature of ferroelectric Ising superlattices have 
been investigated in the context of the mean-field theory [n! and effective-field theory Q, 
and it is found that the dielectric constant has a maximum for a small superlattice period at 
room temperature. It is well known that mean-field theory should not be applied to investi- 
gate fluctuation effects near the phase transition point. Here we adopt effective-field theory 
(based on the Ising spin identities and the differential operator technique) which is superior 
to mean-field theory. It is found that in the low temperature region the larger interfacial cou- 
pling fluctuations, the higher the susceptibility of the superlattice. Remarkably, interfacial 
coupling fluctuations lower the phase transition temperatures when the interfacial coupling 
is stronger than the pseudo-spin coupling in the interior layers. Finally, parameters that im- 
itate BaTiOs/ SrTiOs superlattice are applied to our model. We find that large interfacial 
coupling fluctuations will increase the dielectric constant of BaTiO^/ SrTiO^ structure at 
room temperature, decrease the dielectric constant at high temperature, and further more 
change the critical thickness of the superlattice at which there exists a maximum value of the 
dielectric constant. We postulate that an interfacial disorder such as the interfacial coupling 
fluctuation is one of the reasons for dielectric enhancement in the ferroelectric multilayer. 
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II. MODEL AND FORMULATION 



We consider a ferroelectric superlattice composed of two different components(A and B) 
with an interface (I). L^, L;,, and Lj„ are the thicknesses respectively of A , B, and the 
interface in a unit cell of the superlattice. Each layer is defined on the x-y plane and with 
pseudo-spin sites on a square lattice(see Fig. 1). As in Refs. jll[ 12|, we consider only one 



interfacial layer (Lj„ = 1). The coupling constants between pseudo spins in the interfacial 
layer and in the interior layer can be different from that between spins within the interior 
layers. The system is described by the Ising Hamiltonian with a transverse field, 

H = — JijS^Sj — JinS^S^ — ^iS:[ — 2fiE S^, (1) 

<ij> <mn> i i 

where f2j is the transverse field. S'f, S'f are components of spin-1/2 operator at site i, Jij is 
the coupling constant between the nearest neighbor pseudo spins within the component A 
or B, Jin is the nearest neighbor coupling constant in the interfacial layer and that between 
the interface and the interior layers of component A or B. Jj„ is assumed to be randomly 
distributed according to the independent probability distribution function p{Jin)- /i is the 
effective dipole moment, and E is the applied electric field. The parameters Jij and fli are 
taken as: 

; J. for ,.j e component A 

Jij — \ (2) 

Jfe for i,j & component B, 



for i G component A 



Qi 



^in for i G the interface I (3) 
Qb for i G component B. 
For the ferroelectric material with the first-order phase transition, a four-spin interaction 
term I?! must be included in the Hamiltonian (1). Here, we will focus our attention on 
the effects of the coupling fluctuation in the interface, and the four-spin interaction term 
is not considered for simplicity. The average values of pseudo spins in each layer of the 
superlattice can be derived from the effective-field theory with correlations. For instance. 
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for La — Lb — S, we have 



= « Sf »r=< [cosh{l V Jin) + 2 < 5*8 > sm/i(i V Jin)]- 
[cosh{^ Ja) + 2< Si> sinh{\ y Ja)Y- 

[C0sh{\ V 4) + 2 < -52 > Sinh{\ y -^a)] >r - fix, ^a)\x=0: 

R2 ^ < SI >= [cos/i(| y J„) + 2 < 5"! > sm/i(| y J„)]. 

[cos/i(| y J„) + 2 < ,^2 > sm/i(| y J„)]^. 

[cos/i(| y + 2 < ^3 > sm/i(i y Ja)]-/(a:, fia)U=o, 
R3 = « SI »r=< [cosh{l y Ja) + 2 < 5*2 > sinh{l y J^)]. 

[cosh(i y J„) + 2 < ^3 > sinhil y Ja)]^ 

[cos/i(| y Jin) + 2 < 5*4 > sin/i(i y Ji„)] >r ./(x, Oa)|a;=o, 
R4 = << -S"! >>r=< [cosh{^ y Jin) + 2 < 5'3 > y Jj„)]. 

[cos/i(| y Jin) + 2 < 5'4 > sm/i(| y Jin)]^. 

[cos/i(| y ,/,„) + 2 < 5'5 > sinh{^ y Jin)] >r - fix, Qab)\x=o, 
i?5 = « 5| »,.=< [cos/i(| y Jia) + 2 < 5*4 > sm/i(i y Jin)]. 

[cosh{l y Jb) + 2 < 5*5 > sm/i(| y Jf,)]^. 

[cos/i(| y Jb) + 2 < S'e > y J^)] >r ./(x, 06)U=o, 

i?6 = < -S"! >= [cosh(l y Jb) + 2 < 5'5 > sm/i(| y J^)]. 

[cos/i(| y Jb) + 2 < 5'6 > sm/i(| y J^)]^. 

[cos/i(| y Jb) + 2 < ^7 > sm/i(| y Jb)]./(a;, (]b)U=o, 
i?7 = << >>r=< [cosh{^ y Jb) + 2 < 5*6 > smh{\ y Jb)]. 

[cos/i(| y Jb) + 2 < S'v > sm/i(| y .h)f. 

[cosh{\ y Jin) + 2 < ^8 > sinh{\ y Ji„)] >r ./(x, Ob)]|^=o, 
i?8 = « .51 »r^< [cosh{\ y Jin) + 2 < ,^7 > sinh{\ y Ji„)]. 

[cos/i(| y Jin) + 2 < 5'8 > sm/i(| y Ji„)]^. 

[cosh{\ y Ji„) + 2 < 5"! > sm/i(i y Jin)] >r ■f{x,^ab)\x=0, 



where function f[x) is defined by 



X . ,1 



(4) 



/(a;, ^i) = — tanh(-x + a^). (5) 

To proceed further, we have to approximate the thermal multiple correlations on the right 
side of Eq. (4). We shall use the Zernike decoupling approximation, 

< SfS'...S'kSf >«< S^ >< S' > ... < S^ >< Sf > . (6) 
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In order to describe the coupling fluctuations in a simple way, we assume that the distri- 
bution of Jin is taken to be 

P{jin) = l^^^iJin - Jin " AJj„) + 5{Jin - Jin + AJj„)], (7) 

and the parameter 5j„ (which is introduced to describe the magnitude of the coupling fluc- 
tuation in the interface) is defined as 

5in = (8) 

The symbol <...>,. in Eq. (4) denotes the average over random bonds. These random 
bond averages are given by 

< COsh{-\/Jin) >r = COsh{\/ JinSin)cOsh{\y Jin) , 

[9) 

< sinh{\/Jin) >r = cosh{\/Jin5in)sinh{\/Jin). 

The equations for Ri in (4) ,where i runs over all layers in one period of the superlat- 
tice, form a set of nonlinear simultaneous equations from which each Ri can be calculated 
numerically. The average polarization of the superlattice can then be obtained as 

L 

= ^2/.i?,/L, (10) 

with L = La + Lh + Lin- When the applied electric field E is varied, the average susceptibility 
of the superlattice is obtained numerically from 

dP 

X=g^\E=0. (11) 

By changing the value of 6in, we can investigate the effects of the interfacial coupling 
fluctuations on the susceptibility and the polarization of the ferroelectric superlattice. The 
numerical results and discussions are given in Sec. III. 
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III. RESULTS AND DISCUSSIONS 



Effects of interfacial coupling and the thickness on the ferroelectric multilayer have been 
studied in detail [J, Q]. Here, we focus our attention on the effects of the fluctuations 
of the interfacial pseudo-spin coupling. We first fix the thickness of the superlattice to 
investigate the above effects. In Fig. 2 and Fig. 3, we take La = = 5, Ja = 2Jb, and 
Qa = = O.SJfe, where is taken as the unit of energy. Fig. 2 and Fig. 3 are plotted 
for weak and strong interfacial pseudo-spin couplings respectively. As shown in Fig. 2 



(Jin = V JaJb), for a weak interfacial pseudo-spin coupling, the fluctuations of the interfacial 
coupling result in an increment of the susceptibility of the ferroelectric superlattice only in 
the low temperature region. The phase transition temperature is almost constant as the 
interfacial coupling fluctuation is increased. This is reasonable, because the phase transition 
temperature is mainly determined by the component A which has a stronger pseudo-spin 
coupling than in B and in the interface. The average polarization of the superlattice has 
a slight decrease as the interfacial coupling fluctuation is increased (See Fig. 2(b)). The 
effects of the interfacial coupling fluctuations are more pronounced for a strong interfacial 



coupling {Jin = 3V JaJb)- In Fig- 3(a), the susceptibility of the superlattice will increase 
greatly with increasing 5j„ below the transition temperature, and the peak positions of the 
susceptibility shift to lower temperatures. The dependence of the superlattce polarization 
on the temperature for selected values of are plotted in Fig. 3(b). We can thus see that 
the interfacial coupling fluctuations play an important role in the dielectric properties and 
phase transition temperatures of ferroelectric superlattices for strong interfacial pseudo-spin 
couplings. 

In Fig. 4, we fix the interfacial coupling fluctuation to investigate the effects of the 
superlattice period on the dielectric properties. With the increase of the period of the 
superlattice, it is found that the peak value of the susceptibility decreases and the peak 
position shifts to higher temperatures. And a ferroelectric-ferroelectric phase transition 
occurs for a large period of the superlattice at low temperature, which is also observed by 
use of the mean-field theory in Ref. 

In order to study the interfacial coupling fluctuation effects of real physical systems, we 



consider a BaTiO^/ SrTiOg superlattice with the parameters in Fig. 1 chosen as 11, Q| 
Ja = 26AK, Via = OMK, Jb = 2AK, Vt^ = 87K, J,„ = y/JJ'b, and ili^ = yAUh. We also 
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assume that the interfacial couphng fluctuations are described by Eq. (7). The susceptibihty 
curves are plotted in Fig. 5. When 5j„ = 0, there exists a peak value of the susceptibility 
around = = 4 at room temperature, and our result at Sin = recovers that of 
Ref. 12 . For large fluctuations of the interfacial coupling = 6.0), the peak position of 
the susceptibility occurs at a large value of the period at room temperature. However, as 
shown in Fig. 5(b), the susceptibility at higher temperatures will decease with the increase 

of Sin- 
In summary, the susceptibility (and its peak value) of the superlattice will increase with 
greater interfacial coupling fluctuations in the low temperature region. When the interfa- 
cial coupling is strong, the effects of the interfacial fluctuations on the susceptibility and 
polarization are more pronounced than for weak interfacial couplings. The peak position of 
the susceptibility shifts towards a large period for large interfacial coupling fluctuations at 
room temperature. We conclude that interfacial coupling fluctuations is one of the reasons 
for dielectric enhancement in the ferroelectric multilayer structure. 
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FIG. 1: Schematic illustration of the model of a ferroelectric super lattice. The interfacial layer is 
marked by hollow symbols. 
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FIG. 2: The dependence of the susceptibility(a) and the mean polarization(b) of the ferroelectric 

supcrlattice on the temperature for a given weak interfacial coupling Jj„ and different interfacial 
coupling fluctuations Sin- 
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FIG. 3: The dependence of the susceptibility(a) and the mean polarization(b) of the ferroelectric 
superlattice on the temperature for a given strong interfacial couphng and different interfacial 
couphng fluctuations 6^- 
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FIG. 4: Plot of the susceptibility versus temperature of the ferroelectric superlattice with different 
super lattice periods. 
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